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The local-field (LF) problem of a finite-size dipole emit- ter radiating inside a lossy inhomogeneous structure 
is a long-standing and challenging quantum optical problem, and it now is becoming more important due to 
rapid advances in solid-state fabrication technologies. Here we introduce a simple and accurate quasi-normal 
mode (QNM) technique to solve this problem analyti- cally by separating the scattering problem into contribu¬ 
tions from the QNM and an image dipole. Using a real- cavity model to describe an artificial atom inside a 
lossy and dispersive gold nanorod, we show when the contri- bution of the QNM to LFs will dominate over the 
homo- geneous contribution. We also show how to accurately describe surface scattering for real cavities that 
are close to the metal interface and explore regimes when the surface scattering dominates. Our results offer 
an intui- tive picture of the underlying physics for the LF problem and facilitate the understanding of novel 
photon sources within lossy structures, (c) 2015 Optical Society of America 
OCIS codes: 240.6680, 160.4236, 270.5580. 


The local-field (LF) problem in quantum optics has 
a long, rich history and is now becoming more impor¬ 
tant due to the development of advanced nanofabrication 
technologies for optical structures and nanophotonic de¬ 
vices. As is well known, the local photon density of states 
(LDOS) largely controls light-matter interactions, in¬ 
cluding the spontaneous emission (SE) of a dipole emit¬ 
ter [T] located in an arbitrary scattering environment. 
The LDOS is related to the photon Green function (GF) 
of the medium, through LDOS oc Im[G(r,r' = r;a;)] [5], 
where G(r, r'; w) is the electric-field response at r from a 
dipole emitter at r'. For example, the Im(G) inside a ho¬ 
mogeneous lossless dielectric is simply the vacuum LDOS 
scaled by the refractive index of the dielectric. A plethora 
of novel linear and nonlinear optical effects have been 
predicted and demonstrated in nanophotonic structures 
engineered to have an LDOS that differs substantially 
from that in uniform lossless media [3]. The vast ma¬ 
jority of the experimentally verified LDOS effects have 
been obtained when dipole emitters are located in or in 
proximity to lossless dielectric structures, or in proxim¬ 
ity to lossy (metallic) struc- tures. In all of these cases, 
the Im(G) is well-behaved and there is typically good 
agreement between experimental results of the modified 
LDOS and model calculations using Im(G). Rapid ad¬ 
vances in experimental nanoplasmonics mean that it is 
now possible to embed dipole emitters within inhomoge¬ 
neous metallic microcavities HE]. This raises the chal¬ 
lenging theoretical problem of how to deal with the fact 
that the LDOS is infinity inside lossy, dispersive struc¬ 
tures for which the imaginary part of the permittivity is 
non-zero,i.e., Im[£(a;)] ^ 0. 

Related theoretical work was motivated by dispersive 
lossy materials in which dipole-emitting atoms or ions 


might be interstitially or substitutionally located. Two 
distinct models were developed, both of which isolated 
the dipole “defect” within a cavity that occupied some 
volume in the vacuum between host atoms. The virtual 
cavity (VG) and real cavity (RG) models differed in how 
the LF is treated [7]: in the VG model, the LF is ob¬ 
tained by averaging the macroscopic field inside the cav¬ 
ity, whereas in the RC model, the LF is obtained by 
self-consistently including the scattering from the cavity 
boundary. Relevant experiments, though limited, suggest 
that for ion-doped lossy materials, the RC model works 
well however comparison between experiment and 
model relies on fitting the effective ” cavity radius” of the 
ion, a rather ill-defined physical quantity. A significant 
practical and conceptual advantage of nanofabricated 
structures that surround finite-size nano-dipoles [UlSlIH] 
is that there is typically no ambiguity associated with the 
size and shape of the RC in which the dipole is located, 
and it is manifestly clear that the RC, versus the VC 
model should apply when dealing with the all-important 
LF problem. 

Most works on the RC model to date deal with the LF 
inside a lossy homogeneous material. It has been shown 
the Born approximation could be used to address the ar¬ 
guably more important LF problem in an inhomogeneous 
lossy structure [9] , though it has been confined to spheri¬ 
cal cavity structures. Recently, the LF effect in the strong 
coupling regime for an atom in a spherical structure has 
also been studied [10]. For non-spherical geometries, the 
finite-difference time-domain (FDTD) [TT] algorithm is 
effective for computing the LF; however, the simulations 
are very computational-time expensive and it is difficult 
to understand the underlying physics of the scattering 
problem. Thus one would like a more efficient theoreti- 
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cal approach to describing the SE rate from emitters lo¬ 
cated inside arbitrarily shaped, lossy material structures. 
It would be especially useful to have a formalism that 
facilitates quick exploration of parameter space, e.g., to 
explore the role of RC size and position, on the overall 
device performance. Such models are also relevant to the 
study of spasers [T2l[T3] and quantum plasmonics [14]. 

In this Letter, we introduce a QNM technique devel¬ 
oped for metallic nanoresonator (MNRs) [15UI^ to ad¬ 
dress the LF problem with a RC model wherein the ra¬ 
dius of the RC is directly associated with the volume of 
the emitter, e,g., an artificial atom or quantum dot (QD). 
The MNR problem is particularly challenging from a nu¬ 
merical perspective because of the large complex per¬ 
mittivity, requiring very small (sub-nm) computational 
gridding and tediously long simulation times. However, 
we show that the LF can be solved semi-analytically for 
lossy cavity structures, even for complex-shaped MNRs, 
if one adopts a QNM approach. The QNMs are the dis¬ 
crete modes of an open system, and localized surface 
plasmons (LSPs) may be directly understood as QNMs 
of the MNRs [T5l[T8], defined as the frequency domain 
solutions to the wave equation with open boundary con¬ 
ditions [T9|[20]. The QNM has been widely exploited in 
various field of physics [EIEQ], and a QNM approach 
is particularly useful when the response of the system 
is dominated by one or several discrete resonances. In 
this work, we study the example of a RC inside a gold 
nanorod. To check the accuracy of our theory, we first 
specify a RC embedded at the center of the nanorod, 
and compute the SE enhancement factor using an es¬ 
tablished FDTD technique m and our semi-analytical 
method. We find an excellent agreement between both 
methods for the LF computation. Our approach not only 
offers orders-of-magnitude improvement in the required 
computational time, but helps to identify the underlying 
physics in an intuitive and easy-to-understand way. Us¬ 
ing this semi-analytic technique, we then investigate the 
role of emitter size, and hnd that, for large RC emitters 
such as a QDs, the contribution of the QNM can domi¬ 
nate over the homogeneous contribution of the gold. In 
addition, we also explore the SE rate for RCs placed near 
the MNR surface (inside the lossy structure), and show 
that, ignoring quantum tunnelling effects, the physics 
can be accurately captured using an image dipole tech¬ 
nique. We again explore regimes when the surface scat¬ 
tering can be the dominant mechanism for SE enhance¬ 
ment. 

In the following we consider a 3D gold nanorod as 
our stereotypical arbitrarily-shaped lossy material. Gold 
MNRs can achieve LSP responses in the visible. Other 
geometries including patch antennas have experimen¬ 
tally been probed to achieve enhanced SE of QDs Un- 
Figure [TJa) shows a schematic of the gold rod with a 
length 100 nm, radius = 15 nm, lying along y-axis 
with a background index rib = 1.5 (£& = n^). For the gold 

MNR, we use the Drude model, e(w) = 1 — with 
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Fig. 1. (a) Schematic of a 3D gold nanorod in a homogeneous 
background (n^ = 1.5). (b-c) Spatial profile of the QNM, fc with 
LDc/27r = 324.981 THz (1.344 eV); (b) indicates the LF problem 
with a RC near the center of the nanorod along the y-axis, and (c) 
shows the same RC but 1.5 nm away from the top surface. The 
radius and permittivity of the cavity are rc and Ec, respectively 
(see text). 


parameters: uip = 1.26x10^® rad/s, 7 = 1.41x10^^ rad/s. 

In Ref. [ 22 , Tomas showed that the LF could be de¬ 
scribed by a combination of both the volume averaged 
homogeneous contribution (using the homogeneous GF 
for gold) G}),j.(rd, r^;; w) (r^ is the position of an electric 
dipole at the center of the RC) and the scattering contri¬ 
bution, GfQj,(r(i, r^; w), due to the inhomogeneity of the 
lossy structure; so the LF GF is given by jH] 

Gioc(rd,rd;w) = G!),c(rd, r^; w) -f Gio^(rd, r^; w), (1) 

where G^=^(r(i, r^; w) = ( 2 ^^) G®‘=(rd, r^;w) and 

G^^{rd,rd;u}) is the scattered GF without the RG. Us¬ 
ing a Born-expansion method, Dung et ai, |5] showed 
direct agreement with the Tomas LF formula. The ho¬ 
mogeneous contribution is obtained analytically with the 
two-layer spherical model, yielding a I/r^ dependence at 
a small cavity radius: Im[n-Gj),^-h] « Thus 

one recognizes that the LF problem is essentially solved 
if one can obtain though this is far from trivial, 

except for very simple shapes. 

With regards to MNRs, it has been shown the scat¬ 
tering behaviour can be accurately described in terms 
of the QNMs [16ItT5] . and a mode expansion technique 
for the GF has been developed and confirmed to be very 
accurate plfTTj . The QNMs of the system, f^, have com¬ 
plex eigenfrequencies = cop — iq^, and are normalized 
as follows: [I^[ 2 ^ 
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In practice this normalization is calculated within a 
finite-size computational domain, large enough to cap¬ 
ture the evanescent field contribution of the QNM m- 
An alternative derivation of the normalization scheme is 
presented in Ref. [18] . The corresponding GF with equal 
space arguments within the MNR can then be written 
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as an expansion of the QNMs [MHz], 


G'^(r,r;a;) = ^ 


- w) 




( 3 ) 


For the gold nanorod of interest, we actually only require 
the dominant LSP mode, and thus obtain a single mode 
version of Eq. © as Gj(r,r;a;) = fc(r)fe(r)- 

One problem with the single mode expansion is that 
multi-modal contributions near the surface are not cap¬ 
tured. Such effects becomes important when one is a few 
nm from a lossy structure interface, causing a divergence 
in the LDOS due to quasi-static coupling (Ohmic heat¬ 
ing). However, boundary effects can be well described 
by an image dipole [2l[2l]. Thus, we separate the to¬ 
tal scattered contribution into the contributions of the 
QNM and image dipole: G^^^ = -b Gf„"^F For cavity 
positions sufficiently far from an interface, e.g. 10 nm, 
then one can safely use Gjq^ = for the LF problem. 
For convenience, we define the total enhancement of SE 
(LDOS) rate as Fd = F ^where are the 

QNM, image dipole and homogeneous contributions, re¬ 
spectively, and are given by ^ 

with G®/'^ = and G° = G}),^, respectively (n^ 

is the unit vector of polarization). Here Gq is the free 
space GF, and Im[nd • Go(rd,rd;w) • n^] = uj^/{6ttc^). 

Given sufficient computational resources, it has been 
shown that FDTD works very well when dealing with 
the LFP using a grid size of (^^)^/^ [TT]; the SE rate 
enhancement can be obtained by injecting an electric 
dipole at the center of the RC, and computing the to¬ 
tal response of the dipole. This yields the total LF SE 
emission rate from FDTD, 


pFDTD 


(ra,w) 


Im {tid ■ G^P™(rri,rd;a;) ■ n^} 
Im{nd ■ Go(rd,rd; w) • rid} 


Our goal is to derive Fd analytically, and to first con¬ 
firm its accuracy against a brute-force numerical simu¬ 
lation of Then we use our analytical approach 

to explore some new physics related to LF effects with 
metallic nanorods. 

With a RC, whose radius is much less than the scale 
of the nanorod (for which its effect on the QNM is neg¬ 
ligible), then the QNM contribution is given by [22] 

Gf^^{rd,rd;uj) = {^ ^^e + e ) 


which is obtained analytically from Eq. © with the nor¬ 
malized QNM fc. So numerically, one must first compute 
the QNM, but then, unlike FDTD, Gj)}^ is known for all 
spatial points within the scattering geometry (MNR). 
Furthermore, using a quasistatic approximation-the vari¬ 
ation of electric field produced by the image dipole is 
negligible over the region of the RC, the effect of image 
dipole is given by [21125] 

Gfoc (I’d, Yd] uj) = ^ Gd(rd, -Vd] w), (6) 

Ze + Sc 


with the image dipole GF Gd(r, —rd;u;) = 

“y P“i- 

tion r inside the lossy layer (e), where ± correspond to 
the dipole polarized normal and parallel to the interface, 
respectively; and 1 is the unit dyadic. Below, we will 
concentrate on a y-polarized electric dipole near the top 
surface, so the -I- sign will be used in Eq. and the 
enhanced SE rate from the image dipole is given by 


F = a 
^ \2yduj 


Im 


(sfc - e) 


{sb + £){sc + 2e) 
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where yd is the distance from the center of the cavity to 
the flat surface. A similar image dipole method has been 
used recently [23] to account for the boundary scattering 
effect of the two layer dielectric structure, though we 
find that our method is more accurate as is shown in 
Fig. ©a). 

We now have, in analytic form, the total LF GF for 
the RC inside the lossy MNR, 
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anal 

loc 


= G 


+ G 


qsc 

loc 


+ G 


dsc 
loc ’ 


( 8 ) 


which separates the contributions into (i) a homogeneous 
contribution, (ii) a QNM contribution, and (Hi) a sur¬ 
face contribution. Subsequently, we have an analytical 
prescription for obtaining the total local-field Fy. 
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Fig. 2. LF effect: a y-polarized dipole at 0.5 nm away from the 
center of the nanorod as is shown in Fig. ^b); Tc = 0.62 nm for 
(a-b), and Vc = 3 nm for (d). (a) Full FDTD dipole calculations of 
^FDTD solid), and the analytical result of Fy (green dashed); 
(b) Difference between and Fy (orange solid), and the 

QNM contribution, (c) Ratio of the F^ to the F^ versus cavity 
radius: on resonance (a; = lJc, red solid), redshifted (cd = 1.3 eV, 
orange dot-dashed) , and blueshifted {uj = 1.4 eV, blue dashed), 
(d) Fy (red solid) and the F^ (blue dashed). 


In Fig.©b/c) we show the spatial profile of the dipole 
mode of the LSP for the MNR with resonance around 
1.34 eV [17] . and we are interested in the frequency 
regime around the LSP. The FDTD [2S] calculation of 
the homogeneous contribution is described in Ref. m, 
and it can also be calculated analytically for a spherical 
RG [27l|28]. In Fig. ©a), the homogeneous contribution 
to the enhanced SE rate (i.e., Fy) is shown by the green 
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dashed line for a cavity with radius rc = 0.62 nm and 
refractive index ric = 5 (we use this value for the rest of 
our paper, but it can easily be changed). The FDTD cal¬ 
culation of the total enhancement (red solid) for 

the same cavity located 0.5 nm away from the center of 
the nanorod along the y-axis (Fig. [l](b)) shows a bump 
around the LSP; it is found that the depolarization effect 
of the image dipole, Fy, is negligible at this position. Due 
to the position of the resonance bump in Fy, we expect 
it should be an indication of the LSP. An independent 
calculation of the scattering contribution via the QNM 
using Eq. (O (blue dashed line in Fig.[2jb)) indeed shows 
excellent agreement with the difference between the to¬ 
tal and homogeneous contribution of the enhancement, 

pFDTD _ T^O 

y y ■ 

In Fig. He) we show the ratio of Fy to Fy as a function 
of the radius of the RC. It can be seen as Tc increases, 
the QNM contribution of the LSP begins to dominate 
over the homogeneous contribution. Figure H^) shows 
both Fy and Fy for a RC with radius rc = 3 nm; the 
contribution of the LSP is now larger than the homoge¬ 
neous one, and it also shows a strong asymmetry around 
the LSP resonance. 

The above observation shows that our formulation, 
Eq. (H), works well for computing the LF effect, when 
the RC is at the center of the rod. But for a truly ro¬ 
bust technique, it should also work near a surface. Thus 
in Fig. [3K^)j show the enhanced SE rate for a RC 
inside the nanorod but only 1.5 nm away from the top 
surface along the y-axis (Fig. He)): the green solid line 
is the FDTD result F™™, and the magenta dashed 
is Fy, which is seen to underestimate the total SE en¬ 
hancement. Detailed calculation shows that the QNM 
contribution is now less than 0.1%, and this can be ex¬ 
plained by the mode profile of the QNM, for which the 
nodal lines lay around both ends of the nanorod. The or¬ 
ange dashed line in Fig. Ha-) is calculated with Eq. (H, 
and once again Eq. © agrees well with full FDTD cal¬ 
culation, where now it is important to include the scat¬ 
tering contribution from the image dipole. As shown in 
Fig. Hb), the contribution of the image dipole begins to 
dominate over the homogeneous contribution when the 
separation between the cavity and the top surface comes 
into the sub-nm regime, although quantum tunneling ef¬ 
fects near the surface may well also become important 
here. 

Finally, since it is known that MNRs exhibit a large 
amount of quenching, we have also computed the output 
/3-factor, which is the probability that a single photon 
will be emitter into an output radiation mode, i.e., away 
from the MNR. For a 2 nm RC at the MNR center posi¬ 
tion, the /3-factor is around 9%, which is certainly large 
enough to get light in and out with reasonable efficien¬ 
cies. 

In summary, we have studied the LF problem for a 
finite-size lossy structure. Using a QNM technique, we 
have developed a semi-analytical formulation to compute 
the LF effect for a finite-size emitter inside a MNR. The 
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Fig. 3. LF SE effect for Vc = 0.62 nm. (a) Enhancement at 
1.5 nm away from the top surface as is shown in Fig. He): 
magenta dashed/green solid/orange dashed/cyan circles/ are 
F°/F™™/Fy iva Eq. (O/Fy but using Fy given in [2l], re¬ 
spectively; (b) enhancement at uje versus the distance hy (from 
the top surface to the center of the RC): the red solid/blue dashed 
are Fy and the total calculated via Eq. respectively; the vertical 
dashed line is the regime of (a). 

accuracy of this approach is first confirmed by compar¬ 
ing with an established FDTD technique. The QNM LF 
approach enabled us to identify when one gets a signif¬ 
icant contribution to the LF from the LSP resonance. 
We have also introduced an image dipole technique to 
obtain surface contributions for RCs that are near lossy 
surface. These LF models are of fundamental importance 
in quantum optics and could have applications for LSP 
lasing and spasing. 
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